Quantum disorder in the two-dimensional pyrochlore Heisenberg antiferromagnet 
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We present the results of an exact diagonalization study of the spin- 1/2 Heisenberg antiferromagnet 
on a two-dimensional version of the pyrochlore lattice, also known as the square lattice with crossings 
or the checkerboard lattice. Examining the low energy spectra for systems of up to 24 spins, we find 
that all clusters studied have non-degenerate ground states with total spin zero, and big energy gaps 
to states with higher total spin. We also find a large number of non-magnetic excitations at energies 
within this spin gap. Spin-spin and spin-Peierls correlation functions appear to be short-ranged, 
and we suggest that the ground state is magnetically disordered. 



I. INTRODUCTION 

Geometrically frustrated magnets have attracted the 
interest of both experimentahsts and theorists in re- 
cent years 0, particularly because of their unusual low- 
temperature properties. Fluctuations enhanced by frus- 
tration can destroy long range order, giving rise to dis- 
ordered ground states even in systems with two or three 
space dimensions. Present understanding of such states 
remains incomplete. In this context it is useful to study 
a variety of models, so as to establish the types of be- 
havior that are possible and to isolate good instances of 
each type. 

Examples of geometrically frustrated lattices include 
the triangular and kagome lattices in two dimensions, 
and the pyrochlore lattice in three dimensions. All these 
are based on frustrated units ~ triangles in two dimen- 
sions and tetrahedra in three dimensions - and in the 
last two examples these units are combined in a site- 
sharing arrangement. Properties of the spin-1/2 Heisen- 
berg antiferromagnet with nearest-neighbor interactions 
have been studied for each of these lattices. In the case of 
the triangular lattice, the three sublattice Neel order of 
the classical ground state is believed to be stable against 
quantum fluctuations |^,^ . By contrast, in the case of the 
spin-1/2 kagome Heisenberg antiferromagnet (KHAF), 
analysis of the low-energy spectra for finite clusters in- 
dicates a disordered ground state |^-||. Calculations for 
the three-dimensional spin-1/2 pyrochlore Heisenberg an- 
tiferromagnet present a great challenge; an approximate 
approach |^ yields a disordered ground state in that case, 
too. 

In this paper, we describe properties of the spin-1/2 
Heisenberg antiferromagnet on a two-dimensional ver- 
sion of the pyrochlore lattice, also known as the square 
lattice with crossings or the checkerboard lattice. We 
denote this model by SLWC. The lattice is illustrated 
in Fig.[^. The frustrated plaquettes (those with bond- 
crossings) have the connectivity of tetrahedra, and are 
corner-sharing, as in the three-dimensional pyrochlore 
lattice. The classical Heisenberg model on this two- 




FIG. 1. The two-dimensional pyrochlore lattice, or square 
lattice with crossings. Circles indicate spins of the Heisenberg 
model, and lines denote exchange interactions, all of equal 
strength. 



dimensional lattice is known to have disordered ground 
states , while for the quantum Heisenberg model with 
general spin S it has been proved |ll| that all ground 
states have total spin zero. From ca culations presented 
below of the low-energy spectra for clusters of N sites, 
with N taking even values from = 10 to = 24, we 
conclude that the ground state is disordered at S = 1/2. 

The SLWC joins a rather short list of two-dimensional 
Heisenberg spin systems known to have disordered quan- 
tum ground states. Other examples, apart from the 
KHAF, involve frustration which is not geometric but 
arises either from further neighbor interactions or from 
multiple spin exchange. 

The effects of further-neighbor interactions have been 
studied extensively for the Heisenberg model on the 
square lattice. With nearest-neighbor exchange Ji and 
second neighbor exchange J2, Neel order is destabilized 
by large ground state degeneracy in the classical model 
at an exchange ratio of J1/J2 — 2 |lj]. With spin-1/2, it 
is likely that the ground state is magnetically disordered 



1 



and has a gap for all excitations [p^|']2^] . The SLWC and 
the Ji — J2 model at J1/J2 = 2 are quite closely related, 
in the sense that the Ji — J2 model has second-neighbor 
exchange interactions which, compared to those in the 
SLWC, are twice as numerous but of half the strength. 
It has been noted very recently |^ that a linear spin- 
wave analysis suggests a disordered ground state for the 
SLWC, as in the Ji — J2 model. 

The consequences of frustration by multiple spin ex- 
change have been examined for the spin- 1/2 Heisenberg 
antiferromagnet on the triangular lattice p^|2^. Com- 
petition between two-spin and four-spin exchange leads 
to two disordered spin states: one has an energy gap for 
excitations, while the other is similar to the KHAF in 
having many low-energy excitations. 

II. MODEL AND METHODS 




FIG. 2. The cluster shapes used in the exact diagonaliza- 
tion calculations. 



The Hamiltonian we study is 



Ti. = 2 Si • Sj, 



(2.1) 



where the sum is taken over pairs of spins i and j joined 
by bonds on the lattice of Fig.^ and the Si are spin-1/2 
operators. The magnitude of the exchange interaction in 
this Hamiltonian is the same as that used for the KHAF 
in Ref . , and is a factor of two larger than that used in 
References and 

We may also express the Hamiltonian as a sum over 
tetrahedra, equivalent to square plaquettes with cross- 
ings, by writing 



(2.2) 



where a labels the frustrated plaquettes and = Sq + 
Sfe + Sc -l- Sd is the total spin of plaquette a, to which the 
sites a, b, c and d belong. Expressing the eigenvalues of 
JqP as J(J+1), the quantum number J takes the values 
0,1 or 2. Since each spin is a member of two plaquettes 
with crossings, the operators Jq are not all independent. 
In consequence, it is not possible to construct a state in 
which J = for every a. Nevertheless, we show below 
that J = dominates and that J = 2 has very small 
weight. 

We use symmetry under spin rotations to block- 
diagonalize the Hamiltonian, and employ the Lanczos al- 
gorithm to calculate the low-lying eigenvalues for clusters 
of up to 24 spins. We find eigenvectors by inverse iter- 
ation. We test the accuracy of an eigenvalue obtained 
after j steps of the Lanczos algorithm by performing an 
additional j/10 steps and then recalculating the eigen- 
value: we report only results for which the eigenvalue 
estimate changes by less than 10~^ under this test. 

As a check of the code, we study 12 and 21 site clusters 
for the spin-1/2 triangular lattice Heisenberg antiferro- 
magnet and the KHAF, examined previously in Refer- 
ences [prl and B. We reproduce the energy levels and 



degeneracies to the accuracy reported in each of those 
papers. 

The clusters we study for the SLWC are shown in 
Fig. I For some cluster sizes, several inequivalent shapes 
are possible. In these cases we choose the shape for which 
the ground state energy is lowest, on the basis that these 
shapes presumably have the least frustrating boundary 
conditions. In particular (with the exception of the 14- 
site cluster) the clusters of Fig. || leave unfrustrated the 
two antiferromagnetically ordered states found in the Ji- 
J2 model, at small and large J2/J1 respectively ]l^. 



III. RESULTS 

A. Energy level spectra 

The low-energy spectra for the 16, 18, 20, and 24 site 
clusters are displayed in Figures ^ and ^, labeled by the 
total spin quantum number, S. They are also listed in 
Table |. We notice several distinctive features of these 
results. First, the ground state for each cluster is non- 
degenerate, and has total spin zero as expected from the 
rigorous results of Ref. Second, there is no evidence 
in the spectra of the structure which one would expect if 
the system had Neel order in the thermodynamic limit. 
More specifically, for a Neel-ordered system the quantum 
mechanics of the order parameter gives rise to a set of 
low- lying energy levels [psj-pSl, termed quasi-degenerate 
joint states in recent studies^,^. These are separated in 
energy from other states and have an excitation energy 
that varies with the total spin quantum number S and 
the cluster size N as S{S + l)/N. States of this kind 
are clearly identifiable in the spectra of triangular lattice 
Heisenberg antiferromagnet clusters from the same size 
range Their absence in the SLWC suggests that the 
ground state of the SLWC does not possess Neel order. 
Third, as we discuss below in more detail, the number of 
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FIG. 3. The low energy spectra for the 16 and 18 site clus- 
ters for the SLWC. 
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FIG. 4. The low energy spectra for the 20 and 24 site clus- 
ters for the SLWC. 



states lying within the triplet spin gap is large and in- 
creases with cluster size. Interestingly, both the absence 
of quasi-degenerate states and the presence of low-lying 
singlet excitations are also features of the KHAF 

We comment in passing that the data of Table | include 
some initially unexpected degeneracies for excited states, 
involving factors of 3 and 5. A 3- fold degeneracy on the 
16-site cluster can be understood by recalling that each 
crossed square has the same connectivity as a tetrahe- 
dron, which has four axes of 3-fold rotational symmetry. 
Analogously, cycles of length 3 arise from permutations 
of the tetrahedra in the 16 site cluster. We believe that 
other 3-fold degeneracies appear in a similar way, and 
that 5-fold degeneracies on the 20-site cluster are associ- 
ated with translations. 



B. Dependence on cluster size 

The variation with cluster size of the ground state en- 
ergy per spin is illustrated in Fig.|^. The unit cell of the 
lattice contains two sites, and so it is natural to group 
clusters with even and odd numbers of unit cells sepa- 
rately. Doing this, a common extrapolation to the infi- 
nite system can be found, as shown. The ground state 
energy as well as the triplet and singlet gaps are plotted 
versus 1/A^. (This choice is made simply as a convenient 



TABLE I. Lowest eigenvalues given as an average energy 
per bond, total spin, S, and degeneracies for various sam- 
ple sizes on the square lattice with crossings. A '-' indicates 
that the degeneracies of these levels have not been calculated. 
Degeneracies do not include the 25-1-1 ms degeneracy. 



iV 


2{S, ■ S,) 


S 


deg. 


16 


-0.367045 





1 




-0.351695 





4 




-0.329769 





4 




-0.329403 





12 




-0.325200 





6 




-0.320232 


1 


9 


18 


-0.346423 





1 




-0.337917 





4 




-0.334194 





1 




-0.330624 





4 




-0.329751 





1 




-0.328540 





1 




-0.324801 





8 




-0.323410 





4 




-0.320797 


1 


_ 


20 


-0.359897 





1 




-0.350021 





2 




-0.342692 





2 




-0.336797 





5 




-0.332620 





5 




-0.329264 





5 




-0.328645 





5 




-0.327938 


1 




24 


-0.360776 





1 




-0.355258 





1 




-0.347472 





2 




-0.344195 





1 




-0.344142 





1 




-0.341941 





6 




-0.340329 





1 




-0.338475 





2 




-0.337994 





6 




-0.337811 





6 




-0.336385 





1 




-0.335848 





2 




-0.334265 





6 




-0.333261 





3 




-0.332998 


1 





3 




-1.15 



1/24 1/20 1/16 1/14 

1/N 



FIG. 5. The ground state energy per spin plotted versus 
1/N, where A'' is the number of spins the system. The two 
Unes shown are guides to the eye and are drawn through the 
data for the larger cluster sizes, with an even (diamonds) and 
odd (squares) number of unit cells respectively. 
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FIG. 7. The gap to the first excited state above the ground 
state is plotted versus inverse system size. For all clusters, 
the first excited state is also a singlet state. 



C. Moments of J„ 
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The fact that the Hamiltonian can be expressed as in 
Eq. 2.2, in terms of the total spins of the plaquettes 
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FIG. 6. The gap to the first excited state with higher total 
spin plotted versus inverse system size. 



way of separating points along the horizontal axis.) The 
cluster sizes studied here are small, and there are clear 
difficulties in drawing quantitative conclusions from the 
data. Nevertheless, a definite indication emerges that the 
ground state is magnetically disordered. 

The dependence of the spin gap on cluster size is pre- 
sented in Fig.||. The fact that the spin gap appears to 
remain non-zero in the thermodynamic limit is, of course, 
in contrast to the size-dependence of the longest wave- 
length magnon energy in a system with Neel order. It 
suggests again that the SLWC does not have Neel order. 
A non-zero spin gap is also believed to be a feature of the 
KHAF, but the spin gap for the KHAF is substantially 
smaller than that measured for the SLWC. For example, 
for the 21 site cluster on the kagome lattice, Leung and 
Elser report a spin gap of size 0.5573 Q (scaled to accom- 
modate our factor of 2 in the Hamiltonian), compared to 
a spin gap of size 2.00 for our 24 site cluster. Similarly, 
the Ji — J2 model at J1/J2 = 2 has a non-zero spin gap, 
but one that is smaller by a factor of about 2 at iV = 24 
than for the SLWC H. 



with crossings, leads naturally to the question of what 
component the ground state wavefunction has in sectors 
with each of the possible values for the quantum number 
J for a given . Expanding the ground state wavefunc- 
tion IV') as 



ao|0)+ai|l) + a2|2), 



(3.1) 



where J^jJ) ^ J{J + l)\J) for J = 0, 1 or 2, we de- 
termine the coefficients aj by measuring {'ip\3'^\ip) and 
{ip\3'^\ip). Performing this calculation on the ground 
state for the 24 site cluster, we obtain (J^) = 0.83534 
and (J^) = 2.01968 for every plaquette a. This yields 
|aop = 0.611, |ai|2 = 0.374, and jaap = 0.015. For com- 
parison, a covering of the lattice with nearest-neighbor 
singlets yields |aop = 0.25, jaip = 0.75 and |a2p = 0.0. 
Since our measured value of |a2p is small, we conclude 
that the ground state lies almost entirely within the sub- 
space of states in which every plaquette with crossing 
contains a nearest-neighbor singlet. In addition, it is 
clear that resonance between different such singlet cov- 
erings must make a central contribution to the ground 
state energy. These results may serve as a constraint in 
future attempts to construct an approximate dimer de- 
scription of the ground state, for instance along the lines 
ofRef. M. 



D. Non-magnetic states in the spin gap 

The presence of a large number of non-magnetic ex- 
cited states in the spin gap is a feature that the SLWC 
shares with the KHAF. These states are spread in en- 
ergy across the spin gap, and the lowest of them defines 
the smallest excitation energy in the system: the singlet 
gap. The dependence of the singlet gap on cluster size 



4 



10 



12 



14 



16 1 
N 



8 20 



24 



FIG. 8. The number of non magnetic (S — 0) states in the 
spin gap plotted versus system size. 



TABLE II. Spin-spin correlations for the 24 site cluster. 
Correlations are measured between spins at various separa- 
tions, r, from spin 9 (as labeled in Figures 9, 10, and 11). 
The nearest-neighbor distance is taken to be unity. 
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is shown in Fig.|^. The gap is broadly decreasing with 
increasing cluster size, but it is unclear from these data 
whether it is non-zero or vanishing in the thermodynamic 
limit. In any case, it is very much larger than the singlet 
gap on kagome clusters of similar size: we find for the 
SLWC with iV = 24 a singlet gap of 0.397, to be con- 
trasted with a singlet gap of 3.02 x 10^^ for the N = 27 
kagome cluster [Q. 

The total number of singlet states within the spin gap 
is large and grows with cluster size. In Figure ^ we show 
the number of states, A^r, in the spin gap for each cluster 
size studied. The large, anomalous jump in A^r for the 16 
site cluster is, presumably, a finite size effect particular 
to that sample. A comparison can be made with results 
for the KHAF described in Ref. §. At TV 24, Aat = 
37 for the SLWC and Aa, = 34 for the KHAF, while 
more generally for the KHAF Ajv grows with (even) N 
as 1.14^. 



E. Ground state correlations 




FIG. 9. A coarse-grained pictorial representation of the 
spin Peierls correlation functions measured for the 24 site 
cluster relative to the bond labeled a. Line widths are scaled 
with the magnitude of the correlation function. Positive cor- 
relations are given a solid line while negative correlations are 
represented with a dashed line. Correlations that are less than 
8% of the on-bond correlation are omitted. The values of the 
correlation functions are given in Table III. 



Finally, we report on the measurements of the ground 
state correlation functions for the SLWC. We have calcu- 
lated spin-spin and spin Peierls correlation functions for 
the 24 site cluster. 

As expected both from the absence of quasi-degenerate 
joint states and from the large spin gap, the behavior 
of the spin-spin correlation function suggests that the 
SLWC does not have Neel order. Values are give in Ta- 
ble ||. Correlations fall off rapidly with distance and 
are much smaller in magnitude than in comparison sys- 
tems which are believed to have Neel order (for example, 
the spin-1/2 Heisenberg antiferromagnet on the triangu- 
lar lattice, studied in a 36 site cluster in References |^] 
and [^). Spin-spin correlations on the SLWC are simi- 
lar in magnitude, but slightly larger than, the spin-spin 
correlations on the same sized cluster for the KHAF . 

An alternative to Neel order as a pattern for symme- 
try breaking is spin Peierls or dimer order. To examine 



this possibility, we measure correlations between bond 
operators Da = • Sj, where the bond a joins nearest- 
neighbor spins i and j. We evaluate the connected cor- 
relation function 

Cp{a,b) = {\{DaDt, + Dt,Da)) - {Da){Db). (3.2) 

The results are displayed in Figures |, |l^ and |ll|, and 
listed in Table HI . The behavior of this correlation func- 



tion is clearly influenced by finite cluster size and periodic 
boundary conditions. It is striking nevertheless that the 
correlation length measured along the long side of the 
cluster is no more than one lattice spacing. The magni- 
tudes of spin Peierls correlations are comparable to those 
reported for the KHAF 
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FIG. 10. Spin Peierls correlations measured relative to 
bond 6 and represented as described in the caption of Fig- 
ure 9. 
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FIG. 11. Spin Peierls correlations measured relative to 
bond c and represented as described in the caption of Fig- 
ure 9. 



TABLE 111. Spin Peierls correlations in the 24 site cluster 
of the SLWC. Spin Peierls correlations are measured between 
the bonds labeled a, b, and c, and all other bonds, y, in the 
cluster as pictured in Figures and0. A bond y is repre- 

sented in the table by a pair numbers which label the two sites 
joined by y. Site numbers are as given in the aforementioned 
figures. 
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IV. DISCUSSION 

In summary, the results described in Sec. Ill combine 
to demonstrate that the SLWC has a ground state with 
unbroken spin rotation symmetry. The lack of Neel or- 
der is signaled by the absence of quasi-degenerate joint 
states, by the magnitude and cluster size dependence of 
the spin gap, and most directly by the spin-spin cor- 
relation function itself. It is less clear whether or not 
the translational symmetry of the model is broken in 
the ground state, although the short range of the spin 
Peierls correlations along one direction in Figures ^ - 
[lT| is suggestive of unbroken translational symmetry. If 
translational symmetry is broken, one might expect ei- 
ther columnar dimerization or plaquette resonating va- 
lence bond states, as considered for the Ji — J2 model 
near J1/J2 = 2 Q |l^Jl^,^. Since only a few examples 
are known of two-dimensional Hcisenberg models with 
magnetically disordered ground states, it is a problem of 
some theoretical interest to decide between these alter- 
natives and, especially, to understand the nature of the 
low-lying singlet excitations. 
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